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Abstract. Ground states and dynamical properties of dipolar Bose-Einstein condensate are analyzed based on the 
Gross-Pitaevskii-Poisson system (GPPS) and its dimension reduction models under anisotropic confining potential. We 
begin with the three-dimensional (3D) Gross-Pitaevskii-Poisson system and review its quasi-2D approximate equations 
when the trap is strongly confined in z-direction and quasi- ID approximate equations when the trap is strongly confined 
in x-, {/-directions. In fact, in the quasi-2D equations, a fractional Poisson equation with the operator (— A) 1 / 2 is involved 
which brings significant difficulties into the analysis. Existence and uniqueness as well as nonexistence of the ground state 
under different parameter regimes are established for the quasi-2D and quasi-lD equations. Well-posedness of the Cauchy 
problem for both equations and finite time blowup in 2D are analyzed. Finally, we rigorously prove the convergence and 
linear convergence rate between the solutions of the 3D GPPS and its quasi-2D and quasi-lD approximate equations in 
weak interaction regime. 
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1. Introduction. Quantum degenerate gases have received considerable interests both theoret- 
ically and experimentally, since the first observation of Bose-Einstein condensate (BEC) with dilute 
bosonic gas in f995. The properties of these ultracold dilute quantum gases are determined by the 
short-range, isotropic contact interactions between the particles, which have been studied extensively. 
For those particles with large permanent magnetic or electric dipole moment, dipole-dipole interac- 
tions are non-negligible, and the dipolar interactions are long-range and anisotropic, different from 
contact interactions. Due to these remarkable properties of dipolar interactions, there have been great 
interests to study dipolar BEC in the last decade. In 2005, the first dipolar BEC with 52 Cr atoms was 
successfully realized in experiments at the Stuttgart University |15j . Very recently in 20ff, a dipolar 
BEC with 164 Dy atoms, whose dipole-dipole interaction is much stronger than that of 52 Cr, has been 
performed in experiments at the Stanford University [21] . These success of experiments have renewed 
interests in theoretically studying dipolar BECs. 

In this paper, we will consider the zero temperature mean-field model of dipolar BEC, the three- 
dimensional (3D) Gross-Pitaevskii equation (GPE) with dipolar interaction in dimensionless form 
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where t is time, x = (a;, y) G K 2 and r = (x, z) = (x, y, z) G R 3 are the Cartesian coordinates, 
ip — ip(r,t) is the dimensionless complex- valued wave function, V(r) is a given real- valued trapping 
potential in the experiments, /3 and A are dimensionless constants representing the contact interaction 
and dipolar interaction, respectively, and £/dip(r) is given as 
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with the dipolar axis n = (711,712,723) G R 3 satisfying |n| = \/n 2 + n 2 , + n 3 = 1. Although the 
kernel £/di p is highly singular near the origin, the convolution is well-defined for p G L P (R 3 ) with 
C^dip * P G L P (R 3 ) (p G (1, 00)) 12 . In the context of BEC, the initial data is usually normalized such 
that ||^(-,0)|| i2 = 1. 
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Denote the differential operators d n = n • V and d nn — d n d n , and notice the identity [2] 
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with <5 being the Dirac distribution, we can re-formulate the GPE (jl.ip as the following Gross- 
Pitaevskii-Poisson system (GPPS) [1 [UJ 
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The above GPPS in 3D conserves the mass, or the normalization condition 
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It was proven [2] that when /3 > and — | < A < /3, there exists a unique positive ground state <E> g 
which is defined as the minimizer of the energy functional, i.e. E$r)(Qg) = miri|| $ y , = i £3.0 ($) and 
the Cauchy problem of the GPPS (|1.4|) - (|1.5p is globally well-posed; otherwise there exists no ground 
state and the Cauchy problem is locally well-posed and finite time blow-up may happen under certain 
conditions [2]. 

In many physical experiments of dipolar BECs, the condensates are confined with strong harmonic 
trap in one or two axes directions, resulting in a pancake- or cigar-shaped dipolar BEC, respectively. 
Mathematically speaking, this corresponds to the anisotropic potentials V(r) of the form: 

Case I (pancake-shaped) , potential is strongly confined in the vertical z direction with 



V(r) = V 2 (x) + r=(x,z)e 



(1.8) 



Case II (cigar-shaped), potential is strongly confined in the horizontal x = (x,y) £ R 2 plane with 

V(r) = V 1 (z) + ^f ) r=(x,z)eM 3 , (1.9) 

where < e <C 1 is a small parameter describing the strength of confinement. In such cases, the above 
GPPS in 3D can be formally reduced to 2D and ID, respectively [TT] . 

In Case I, when e — > + , evolution of the solution ip(r, t) of GPPS (|1.4j) - (ll.5[) in z-direction would 
essentially occur in the ground state mode of L z :— —\d zz + -^x, which is spanned by w e (z) = 
E -V2 7r -i/4 e -5?r [IT1|. By taking the ansatz 
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the 3D GPPS (fl~4|) - ([L5)) will be formally reduced to a quasi-2D equation I [TT] : 
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where x = (x, y), n ± = (ni,n 2 ), d n± = n ± • V, d n±n± = d n± (d n± ), A = d xx + d yy and 
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In addition, as e — > + , tp 2D can be approximated by ip 2 ^ [TT] as : 

1 
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which can be re-written as a fractional Poisson equation [TT] 
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Thus an alternative quasi-2D equation II can be obtained as [11) : 
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Similarly, in Case II, evolution of the solution ip(x, z, t) of GPPS (|1.4[l - (|1.5jl in x = (x, y)-directions 

1 2 , 2 

would essentially occur in the ground state mode of i x := — ^(d xx + d yy ) + , which is spanned 

| x |2 

by w e (x) = e~ 1 TT~ 1 l 2 e~^ T [TTJ[4]. Again, by taking the ansatz 



(x,z)eR J , £>0, 

the 3D GPPS (fl~4)) - (ll~5)) will be formally reduced to a quasi-lD equation [UJ: 
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The above effective lower dimensional models in 2D and ID are very useful in the study of dipolar 
BEC since they are much easier and cheaper to be simulated in practical computation. In fact, for the 
GPE without the dipolar term, i.e. A = 0, there have been extensive studies on this subject. For formal 
analysis and numerical simulation, the convergence rate of such dimension reduction was investigated 
numerically in [51 [S] and a nonlinear Schrodinger equation with polynomial nonlinearity in reduced 
dimensions was proposed in [53]. For rigorous analysis, convergence of the dimension reduction under 
anisotropic confinement has been proven in the weak interaction regime [H][E]> i-e. P — 0(e) in 2D 
and j3 = 0(e 2 ) in ID. However, with the dipolar term, i.e. A ^ 0, there were few works towards 
the mathematical analysis for this dimension reduction except some preliminary results in [T2] where 
different scalings and formaulation were adapted. In fact, our quasi-2D models (II. lip and (|1.15l) and 
quais-lD model (I1.17[) are much easier to be used in mathematical analysis and practical numerical 
computation. 

The main aim of this paper is to establish existence and uniqueness of the ground states and 
well-posedness of the Cauchy problems associated to the quasi-2D equations I and II and quasi-lD 
equation, and to analyze the convergence and convergence rate of the dimension reduction from 3D to 
2D and ID. In order to do so, without loss of generality, we assume the potential Vd(r]) > for n £ R d 
(d = 1, 2, 3). It is natural to consider the energy space in d-dimensions (d = 1, 2, 3) defined as 
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and the unit sphere of X d defined as 

S d = X d f){u £ L 2 (R d )\\\u\\ L 2 { 
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This paper is organized as follows. In Sections 2, 3 and 4, we study quasi-2D equation I (jl.lip . 
II (|1.15p and quasi- ID equation (|1.17[) . respectively. In each section, we first establish existence and 
uniqueness as well as nonexistence of the ground state under different parameter regimes, and then 
study the well-posedness of the corresponding Cauchy problem. In Section 5, we rigourously prove 
the validity of dimension reduction from 3D GPPS (|1.4[) - (|1.5[) to 2D and ID in the weak interaction 
regimes. Our approach is based on a-priori estimates from the energy and mass conservation together 
with the Strichartz estimates. 

Throughout the paper, we adopt the standard notation of Sobolev space and use ||/||£ := J Rd |/(?7)| p d-q 
for p G (0, oo) when there is no confusion about the space M. d , denote C as a generic constant which 
is independent of e, let X* as the dual space of X, and adopt the Fourier transform of a function 
f(rj) G L\R d ) as 

/(£) = / f(v)e- %( - n d V , £eM d . (1.19) 



2. Results for the quasi-2D equation I. In this section, we prove existence and uniqueness 
as well as nonexistence of ground states for the quasi-2D equation I under different parameter regimes 
and local (global) existence for the Cauchy problem. For considering the ground state in 2D, let Cf, 
be the best constant from the Gagliardo-Nirenberg inequality [29], i.e. 

C h := inf "^"^i'" 7 " 2 ^ . (2.1) 

2.1. Existence and uniqueness of ground state. Associated to the quasi-2D equation I (|l.llj) - 
(I1.12p . the energy is 



where 



i|V$| 2 + F 2 (x)|$| 2 + g X t? nlX \*t - ^1*1 V D 



cbc, $ G X 2 , (2.2) 



<?»> = (d n±n± - n 2 3 A) <p 2D , V 2D = U 2D * l*| 2 - (2-3) 
The ground state $ ff G S 2 of (jl.lip is the minimizer of the nonconvex minimization problem: 

Find $„ G S 2 , such that E 2D (<& g ) = min E 2D ($). (2.4) 

For the ground state, we have the following results: 

Theorem 2.1. (Existence and uniqueness of ground state) Assume < V 2 (x.) G L^ C (R 2 ) and 
lim V 2 (x.) = oo, then we have 

|x|— too 

(i) There exists a ground state $ g G S 2 of the system fl.llty - lfl.lSI) if one of the following conditions 
holds 

(Al) A > and - A > - V2^C b e; 

(A2) A < and [3+ ±(1 + 3|2n§ - 1|)A > -y/2irC b e. 

(ii) The positive ground state \$> g \ is unique under one of the following conditions: 
(Al') A > and (3 - A > 0; 

(AS!) A < and (3 + ±(1 + 3|2n 2 - 1|)A > 0. 
Moreover, any ground state is of the form $ 9 = e j8o |$ g | for some constant 9 G R. 

(in) If (3 + ijA(l — 3n|) < — \f2~jrCbS, there exists no ground state of the equation 111, lip . 

In order to prove this theorem, we first study the property of the nonlocal term. 

Lemma 2.2. (Kernel U 2 e D in H~W\) ) For any real function /(x) in the Schwartz space S(R 2 ), we 
have 

U? D * /(£) = f(0 U^(0 = ™l / j^—zds, f G S(R 2 ). (2.5) 
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Moreover, define the operator 



T aa '(f) = d aa >(U e */), a,a' = x,y, 



then we have 



\\T aa >fh < 



V2 
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ll/lb, \\TaafSh < llWIb, 



hence T aa r can be extended to a bounded linear operator from L 2 (R 2 ) to L 2 (M 2 ). 
Proof. From (|1.12|) . we have 



C4 2D (x) 



2V2tt 3 / 2 Jr ^|x| 2 +e 2 s : 



: dS 
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2tt|x 
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This immediately implies that U 2D *g is well-defined for any g € L 1 (M 2 ) p| L 2 (R 2 ) since the right hand 
side in the above inequality is the singular kernel of Riesz potential. Re- write t/| D (x) as [TT] 



using the Plancherel formula, we get 



w' 2 (z)w 2 (z / ) 
v/|x| 2 + (z - z'f 



dzdz' , 



x e 



+ g + 



y R ici 2 + s ; 



rds, 
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which immediately implies f|2 . 5|) . Here c denotes the complex conjugate of c. Concerning T QQ ,<, we only 
need to prove the results for T xx since others are similar. Applying the Fourier transform, we have 



T xx f(0 



/(0 



2,2 

Si 



IsT 
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< 



/(0 



e -^ 2 / 2 ds 
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/(0 
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Thus we can get the first inequality in 
from 



and know that T xx : L 



L? is bounded. Moreover, 



T xx f(0 



/(o 



-e 2 s 2 /2^2 



* J* lsT + * 2 
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1/(01 iai 2 



;ds < ICI 1/(01, 
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we obtain the second inequality in (|2.6j) and know that T xx : H 1 — > L 2 is bounded too. □ 

Remark 2.1. In fact, T aa i is bounded from LP —} LP , i.e., there exists C p > independent of e, 
such that 

\\Taa'{f)\\p < ^ll/llp, PG(1,«>). (2.11) 

T/iis caw &e obtained by Minkowski inequality and LP estimates for Poisson equation. 
Lemma 2.3. For the energy E2d(~) in h2.2\) . we have 

(i) For any $ £ S2, denote p(x) = |$(x)| 2 , then we have 

E 2D m > E 2D {\<&\) = E 2D (Vp) , V$eS 2 , (2.12) 

so i/ie ground state Q g of \2. 2\) is of the form e i6a \^ g \ for some constant 9q £ R. 

(ii) Under the condition (Al) or (A2) in Theorem \2.1{ E 2 D{\fp) is bounded below. 
(Hi) Under the condition (Al' ) or (AS!) in Theorem ] 2. 11 E2d(^/P) is strictly convex. 
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Proof, (i) For any $ £ S 2 , then |$| S 52, and a simple calculation shows 

E 2D (<s>) - E 2D m) = i||V$||l - i||V|$| ||i > 0, $ e S», 



where the equality holds iff 



which is equivalent to 



|V$(x)| = V|*(x)|, a.e. xe 



$(x) = e ie °|$(x)|, for some 8q G K. 

Then the conclusion follows. 

(ii) For yfp = <f> e 5*2, we split the energy E 2 d into two parts, i.e. 



E 2D {$) = + £ 2 ($) = E^y/p) + E 2 (J-p), 



where 



e 2 (Vp) = 



-|V^| 2 + V2(x)p 



/3-A + 3n 2 A, l2 3A ~ 
2 v / 2^£ 4 ^ 



dx, 



with 

Applying the Plancherel formula and Lemma [2T2l there holds 



^( X )p(x)dx = i -2 / ^(OpCO* 



-1 

4-7T 3 



lei 2 + s 2 

Recalling the Cauchy inequality and n\ + n 2 + n 2 , = 1, we have 

- nj\e < (mfc + n 2 6) 2 ~ n 2 |^| 2 < (1 - 2n 2 )|£| 5 



|p| 2 ds^. 



Denoting Co = max{|n||, |1 — 2n 2 |}, we can derive that 



ip 2D (x)p(x) dx 



- 4^3 



/3lpl^=^l H ll. 



/7T £ 



(2.13) 

(2.14) 
(2.15) 

(2.16) 

(2.17) 
(2.18) 

(2.19) 



Hence, E 2 (^fp) can be bounded below by ||p||a- In fact, under the condition (Al), L 



P — A > —\J2iiCb £, we have 



/3 - A + 3n§A 



2V2tt£ 



IIpII 



3\/2n§A 



IIpIII > ~\\p\\ 



2- 



Similarly, under the condition (A2), if A < and n 2 > |, then 

- X -I- 3 

£ 2 (v^) > 



/3- A + 3n 2 A 2 . C 6 2 

IIHI 2 > — dMIa; 



2V2tt£ 



(2.20) 

(2.21) 

(2.22) 
A > and 

(2.23) 
(2.24) 
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and if A < and n| < 5, then 

jp , /3-A + 3n§A „ ||2 3V2(l-2 ra |)A „ ||2 ^ GV 2 

£MVp) > — ^-7s= \\P 2 + 7-7=- p 2 > — H p U- ( 2 - 25 ) 

Recalling the choice of the best constant Cf,, under either condition (Al) or (A2), the energy 

e 2D (Vp) = ExWft + e 2 (Vp) > \wVp\\1 - y IIpIII > o. (2.26) 

(hi) Again, we split the energy as (|2.16[) . It is well known that Ei(<Jp) is strictly convex in p 20 . 
It remains to show that E 2 (yfp) is convex in p. For any real function u E i 1 (R 2 ) n L 2 (R 2 ), let 



H{u) 



+ 3%A , u|2 _ 3A u _ ^ {U , D # u) 



2V2Ve 4 



dx. (2.27) 



Then E 2 (yfp) — H{p). It suffices to show that (p) is convex in p. For this purpose, let ^fp{ = $i £ S 2 
and = $ 2 e 5 2 , for any G [0, 1], consider p g = 6*pi + (1 — 9)p 2 and y/p^ € 5*2, then we compute 
directly and get 

0H{px) + (1 - 9)H( P2 ) - H{ Pe ) = 6(1 - 6)H( Pl - p 2 ). (2.28) 

Similar as (|2.20l) . looking at the Fourier domain, we can obtain the lower bounds for H(pi — p 2 ) under 
the condition (Al') or (A2'), while replacing Cb with in the above proof of (ii), i.e., 

H( Pl -p 2 )>0. (2.29) 

This shows that H(p), i.e. E 2 (y/~p), is convex in p. Thus E2d{-J~P) is strictly convex in p. □ 

Proof of Theorem 12.11 (i) We first prove the existence results. Lemma T2.3I ensures that there 

exists a minimizing sequence of nonnegative function {$™}^ C S 2 , such that lim E 2 D(& n ) = 

n— >oo 

inf £'($). Then, under condition (Al) or (A2), there exists a constant C such that 
4>eS 2 



|V$™|| 2 + ||$"|| 4 + / T/ 2 (x)|$"(x)| 2 dx< C, n> 
Jr 2 



(2.30) 



Therefore $" belongs to a weakly compact set in £ 4 (JR 2 ), iJ^R 2 ), and £y- 2 (R 2 ) with a weighted L - 
norm given by \\<Z>\\l V2 = [J R2 |$(x)| 2 t/ 2 (x)dx] 1 /2. Thus, there exists a $°° e := ff^R 2 ) p| L 2 ^ (R 2 ) p| i 4 
and a subsequence of {<&™}$£L (which we denote as the original sequence for simplicity), such that 

inW, V$"^V$°°, ini 2 . (2.31) 

The confining condition lim F 2 (x) = oo will give that \\$°°\\2 = 1 QHHQ]. Hence $°° <E 5 2 and 

|x|— >oo 

$" — > $°° in i 2 (R 2 ) due to the i 2 -norm convergence and weak convergence of {$ n }^° =0 . By the lower 
semi-continuity of the H 1 - and L^-norm, for E\ in (|2.17p . we know 

< lim inf £d($"). (2.32) 

n— ^cxd 

By the Sobolev inequality, there exists C(jp) > depending on p > 2, such that ||$ n || p < C(p)(\\ V$"|| 2 + 
ll^"^) < C(p)(l + C), uniformly for n > 0. Applying the Holder's inequality, we have 

||($«)2 _ ($°°) 2 || 2 < Ci(||$ n ||i + ||$ n ||i)||$" - $°°|| 2 , (2.33) 
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which shows p n = ($™) 2 — > p°° = (<J>°°) 2 in L 2 (R 2 ). Using the Fourier transform of U 2D in Lemma 
12.21 and (|2.22l) . it is easy to derive the convergence for E 2 in (|2.18p . i.e. 

E 2 (<f>°°) = lim E 2 (<f> n ). (2.34) 

n— too 

Hence 



E 2D {<S>°°) = E^ 00 ) + £ 2 ($°°) < liminf£J 2D ($ n ). (2.35) 

n— >oo 

Now, we see that is indeed a minimizer. For the uniqueness part, it is straightforward by the strict 
convexity of i?2_o(v / P) shown in Lemma 12.31 

(ii) Since the nonlinear term in the equation behaviors as a cubic nonlinearity, it is natural to 
consider the following. Let $ G 52 be a real function that attains the best constant C& [29], then 
$(x) is radially symmetric. Choose $^(x) = (5 _1 $(<5 _1 x) with S > 0, then $5 £ S 2 . Denote ips = 
(3 n±ni — n 2 A±^j (U 2D * \<&s\ 2 ), by the same computation as in Lemma \2. 31 there holds 



47T 3 7 R 3 |£| 2 + S 2 

±_ f (nxa+n 2 6) 2 -n|lgl 2 



|$| 2 (0 



dsd£ 

2 

dsd£. 



4tt 3 (5 2 7 r3 ICI 2 + <S 2 s 2 

Using the fact that $(x) is radially symmetric, |<&| 2 (£) is also radially symmetric, then we obtain 

^l^l 2 rfx = - ^ + ^- 2n | +0(1) p||t as f-O* (2.36) 

Hence, as S —> + , we get 

^(•*) - i (Will + ^ + + + £ ^ x )|$| 2 (x ), x . 

Recalling that ||V$|| 2 = C b ||$|||, we know lim E^d^s) = -00 if /3 + §A(1 - 3n 2 ) < -y/2nC b e, i.e. 

there is no ground state in this case. □ 

2.2. Well-posedness for the Cauchy problem. Here, we study the well-posedness of the 
Cauchy problem corresponding to the quasi-2D equation I (|1.11[) - (I1.12[) . Using the Fourier transform 
of the kernel U 2D in Lemma 12.21 it is straightforward to see that the nonlinear term introduced by 
U 2D behaves like cubic term. Thus, those methods for classic cubic nonlinear Schrodinger equation 
would apply [131 [2^1 [37] . In particular, we have the following theorem concerning the Cauchy problem 
of (fl3B - (Q2l) . 

Theorem 2.4. (Well-posedness of Cauchy problem) Suppose the real-valued trap potential satisfies 
V 2 (x) > for x£l 2 and 

U 2 (x) € C°°(IR 2 ) and L> k U 2 (x) e L°°(R 2 ), for all k e N 2 , with |k| > 2, (2.37) 

then we have 

(i) For any initial data <f>(x.,t = 0) = </>o(x) G X 2 , there exists a Tmax G (0, +00] such that the 
problem \1.11)) - Kl.l2\) has a unique maximal solution <fi e C ([0, Tmax), X 2 ). It is maximal in the sense 
that ?/T m ax < 00, then \\(f>(-, t)\\x 2 — ► 00 when t —} T max . 

(ii) As long as the solution c£(x, i) remains in the energy space X 2 , the L 2 -norm ||<^(-,i)|| 2 and 
energy E 2 r>{(j>(-,t)) in H2.2\) are conserved for t £ [0,T max ). 

(Hi) Under either condition (Al) or (A2) in Theorem \2.1\ with constant C'b being replaced by 
Cb/\\(f>o\\ 2 , the solution of hl.ll}) - fl.l2}) is global in time, i.e., T max = 00. 
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Proof. The proof is standard. We shall use the known results for semi-linear Schrodinger equation 
[T3] . For tfi G X-2, denote p = |</>| 2 and consider the following 

G(<j>, $) := G(p) = \( |0| 2 (d n±n± - n 2 A) (U 2D * |0| 2 ) dx, 



2 

gtf) = S -^Ml = cf, (d nxnx - nlA) (U? D * |0| 2 ), G X 2 
Then the equations (|l.lll) - (|1.12l) read 



id t (j) = 



-~a + f 2 ( x ) 



(f> + /3 \<l>\ 2 (f>-3Xg((f>), xel 2 , *>0, (2.38) 



V2n e 

Sobolev inequality, for ||u|| v 2 + ||«||x 2 < M, it is easy to prove the following 



where /3q = ^ 2 " 3 . Using the L p boundedness of Tj^ (cf. Lemma [2.21 and Remark 12. ip and the 



\\g(u)-g(v)\\ 4/3 <C(M)\\u-v\U. (2.39) 

In view of the standard Theorems 9.2.1, 4.12.1 and 5.7.1 in JTS] and [23 f° r the well-posedness of 
the nonlinear Schrodinger equation, we can obtain the results (i) and (ii) immediately. The global 
existence (hi) comes from the uniform bound for ||0( - ,t)||x 2 which can be derived from energy and 
L 2 -norm conservation. □ 

When the initial data is small, there also exists global solutions [T3JII2]- Otherwise, blow-up may 
happen in finite time, and we have the following results. 

Theorem 2.5. (Finite time blow-up) For any initial data (f>(x,t = 0) = 4>o(x) G X 2 with 
L 2 |x| 2 |</>o(x)| 2 dx < oo, if conditions (Al) and (A2) with constant Cf, being replaced by Gfe/H^oH! 
are not satisfied and assume V^(x) satisfies 2V<2(x) +x- W 2 (x) > 0, and let <j> := (f>(x,t) be the solution 
of the problem U. 11)) . there exists finite time blow-up, i.e., T max < 00, if A = 0, or A > and n| > ^, 
and one of the following holds: 

(i) E2D(<M<0; 

(ii) E 2D (<p Q ) = and Im (/ M2 <j) (x) (x • V^o(x)) dx) < 0; 

(Hi) E 2 d(M > and Im (/ R2 <f> (x) (x • V<?f>o(x)) dx) < - ^/2E 2D (4>q)\\ xcf> \\ 2; 
where Im(/) denotes the imaginary part of f. 
Proof. Define the variance 



a v (t):=a v (<t>(;t))= |x| 2 |0(x, t)\ 2 dx = a x (t) + a y (t), t>0, (2.40) 

where 

a a (t) := a a ((/>(;t)) = I a 2 \^(x,t)\ 2 dx, a = x,y. (2.41) 

For a — x, or y, differentiating (|2.4ip with respect to t, integrating by parts, we get 
d f 

—a a (t)=-i [a4>{x,t)d Q (j){x,t) - a<j)(x,t)d Q 4>{x,t)] dx, t>0. (2.42) 
at J R 2 

Similarly, we have 

i<r«(i) - I [2\d a <t>\ 2 + /?o|0| 4 + 3X\<l>\ 2 ad a {d n±a± - n 2 A)<^ - 2a\0\ 2 d a V 2 (x)] dx, (2.43) 
at Jul 
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where /3 = P ~^ n 1 , V = U 2D * |0| 2 . Writing p = |0| 2 , <p = (d n±n± - n 2 3 A)<p, = (m& + 



— "-lICI an( i noticing that p is a real function, by the Plancherel formula, we have 

(x • v*) rfx = ^ | r2 ^(o v • (^) = ^ £ p(o v • (en^) « 

1 
1 

4^2 



£5 / ( |p| 2 V(£n^) + r^[/ 2 ^ • -Vlpl 2 ) d£ 



|^|Vdx+ — 



4tt3 y R3 del 2 + s 2 ) 2 



dsdt;. 



Denote 



1(f) : =J(0(.,t)) = _ 



using n 5 e [— ™ 2 |£| , (1 - 2n 2 )|£| 2 ], we obtain 



s 2 e- £ s / 2 \p\ 2 



-V2l 



m 2 + s 2 



\/2(l - 2n|) 



dsd£, i > 0, 



V7T£ 



f > 0. 



(2.44) 



(2.45) 



If A = 0, or A > and ri3 > \, noticing XI (t) < in these cases, summing (|2.43p for a = x, y, and 
using the energy conservation, we have 



dt 2 



av(t) = 2 



5- 



|V0| 2 + /3 |^| 4 + g A|0| 2 (x • Vp) - |0| 2 x • W 2 (x) 



<2x 



Thus, 



= 4E 2D {cj>(;t))+3XI(t)-2 [ |0| 2 (2^ 2 (x) +x- W 2 (x))dx 
<m D {<f>(;t)) = 4E 2D (cb ), f>0. 

a v (t) < 2E 2D {(j) )t 2 + a' v (0)t + cry(0), t > 0, 



and the conclusion follows in the same manner as those in 
equation. □ 



for the standard nonlinear Schrodinger 



3. Results for the quasi-2D equation II. In this section, we investigate the existence, unique- 
ness as well as nonexistence of ground state of the quasi-2D equation II (|1.15[) and the well-posedness 
of the corresponding Cauchy problem. 

3.1. Existence and uniqueness of ground state. Associated to the quasi-2D equation II 
(|1.15j) . the energy is 



1 



V*| 2 + V 2 (x)|$| s 



/3-A + 3n|A 4 _ 3A 



2v / 2tt£ 



m 2 v 



dx, $el 2 



(3.1) 



where 



rtx) = (d n±n± -n 2 A)((-A)-V 2 \$\ 2 ). 



(3.2) 
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The ground state $ g £ S 2 of the equation (|1 . 15|) is defined as the minimizer of the nonconvex mini- 
mization problem: 

Find $ g e S 2 , such that E 2 D{^ g ) = min E 2D (<&). (3.3) 
For the above ground state, we have the following results. 

Theorem 3.1. (Existence and uniqueness of ground state) Assume < V 2 (x) € Lf^ c (M. 2 ) and 
lim V 2 (x) = 00, then we have 

|x|— >oo 

(i) There exists a ground state <& g € S 2 of the equation il.15]) if one of the following conditions 
holds 

(Bl) A = and /3 > ~^C b e; 
(B2) A > 7 n 3 = and P~ \> -V2^C b e; 
(B3) A < 0, n\ > \ and [3 - (1 - 3n§)A > -^C h e. 
(it) The positive ground state \<fr g \ is unique under one of the following conditions 
(Bl') A = and B> 0; 
(BS!) A > 0, n 3 = and 8 > A; 
(BSf) A < 0,. n\ > \ and 8 - (1 - 3n 2 )A > 0. 
Moreover, any ground state $ g = e* e ' |<l > g | for some constant 8q £ K. 

(Hi) There exists no ground state of the equation U.15]) if one of the following conditions holds 
(Bl") A > and n 3 ^ 0; 
(BS!') A < andnl <\; 
(B3") A = and 8 < ~V2^C b e. 
Again, in order to prove this theorem, we first analyze the nonlocal part in the equation (|1.15[) . 
In fact, following the standard proof in we can get 

Lemma 3.2. (Property of fractional Poisson equation ll.lSfy ) Assume /(x) is a real valued function 
good enough, for the fractional Poisson equation 

(-A)~ 1/,2 </j(x) = /(x), xel 2 , lim p(x)=0, 

|x|— }OQ 

we have 

p(x) = / 9 f (x ° n dx! = f ^-t) * /, x e K 2 , 

V 2tt|x - x'| V 27r l x l/ 
and the Hardy-Littlewood-Sobolev inequality implies 

IMIp. < C P ||/|U P* = ^- Pe(l,2). (3.4) 
Moreover, the first order derivatives of ip are the Riesz transforms of f and satisfy 

\\d a (f\\ q < Cq\\f\\q, q G (1, oo), a = x, y, (3.5) 
and the second order derivatives satisfy 

\\d aa ,<p\\ q = \\d a ((-A)- 1 ^d a ,f)\\ q <C q \\d a ,f\\ q , ?€(l,oo), a,o/ = x,y. (3.6) 

Remark 3.1. Similar results hold for T aa i defined in Lemma \2.2[ i.e. 

\\T aa .-f\\ v < CpHV/Hp, /or p €(1, oo). (3.7) 
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Since the fractional Poisson operator (— A) -1 / 2 is taken as an approximation of U 2D (|1.12j) . we 
consider the convergence regarding with the derivatives. 

Lemma 3.3. For any real-valued function f G L P (M. 2 ), let 



T e a {f) = d a (U 2D */), R a (f) = ^(-A)- 1 / 2 /, a = x,y, 



(3.8) 



then is bounded from LP to LP for 1 < p < oo with the bounds independent of e. Specially, for any 
fixed f G L P (M. 2 ) with p G (1, oo), we have 



lim . \\T*(f)-R a (f)\\ P = a, pe (l,oo). 



(3.9) 



Proof. We can write R a and as 

R a (f)= K <**f, 
where R a is the Riesz transform and 

1 



T*(f) = K*f, 



^l(x) 



V/2 



2tt|x| 3 ' " a ""' 2V2tt 3 /2 7 r (|x| 2 + e 2 s 2 ) 3 / 2 
It is easy to check that K e a satisfies 

|^(x)| < B|x|- 2 , |V^(x)| < B|x|- 3 



ds, x G 



a = x,y. 



(3.10) 



(3.11) 



Ixl > 0, 



ifj(x)dx = 0, < i?i < i? 2 < oo, 

fli<|x|<_R 2 

for some e-independent constant B. Then the standard theorem on singular integrals |25) implies that 
is well defined for LP functions and is bounded from LP to LP with e-independent bound. 
Thus, we only need to prove the convergence in L 2 , other cases can be derived by an approximation 
argument and interpolation. For the L 2 convergence, looking at the Fourier domain, we find that 



\\T*(f)-RM)\\l = 



< 



1 

I^ 2 
1 

i^ 2 







2 


/ 

/R 2 










2 


/ 


no 




/R 2 







a 

W\ 

1 

7T 



a 

7T 
(1 



-£ 2 S 2 /2 



-i 2 



>iei 



ds 



Notice that for fixed ^ £ G R 2 , the dominated convergence theorem suggests that 



lim — 

r^0+ 7T 



(1-e 



-e 2 s 2 /2) 



lei 2 



= 0, 



(3.12) 



hence, the conclusion in L 2 case is obvious by using the dominated convergence theorem again. Using 
approximation and noticing that L 2 n L q is dense in L p for q G (1, oo), noticing the uniform bound on 
T^ : L p L p for p G (1, oo), we can complete the proof. □ 

Lemma 3.4. For the energy E2d{-) in h3.1\) . the following properties hold 

(i) For any G S2, denote /?(x) = |$(x)| 2 , then we have 



E 2 d($) > E 2D (\<f>\) = E 2D (Vp) , V$e5 2 , 



(3.13) 



so the ground state $> g of h3. 1)) is of the form e ie °\^ g \ for some constant 80 G M. 

(ii) Lf condition (Bl) or (B2) or (B3) in Theorem ] 3. 1\ holds, then E 2 d * s bounded below. 

(Hi) Lf condition (BL) or (BSl ) or (B3 ) in Theorem \3.1\ holds, then E 2 d{^/p) is strictly convex. 
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Proof, (i) It is similar to the case of Lemma 12.31 

(ii) Similar as Lemma [2~B1 for $ € S%, denote p = |$| 2 , we only need to consider the lower bound 
of the following functional 

H(p) = —X f p (d n±a± - n 2 3 A) [(-A)" 1 /^] dx . (3.14) 

JR 2 

Using the Plancherel formula and Cauchy inequality, for A < and n| > 5, we have 

h(p) = ± I \m\ 2 > A f (1 _ 2ni m \ m f # > o. (3 .i5) 

For A > and 713 = 0, it is easy to see H(p) > 0. Hence, assertion (ii) is proven. 

(iii) Similar as Lemma 12.31 it is sufficient to prove the convexity of H(p) in p. For ^fp~\ € S 2 , 
^fp~2 & S 2 and any 9 € [0, 1], denote p e = 0pi + (1 — 9)p 2l we have 

9H( Pl ) + (1 - 0)#(p 2 ) - H(p e ) = 9(1 - 9)H{p x - p 2 ), (3.16) 
where the RHS is nonnegative under the given condition, i.e., H{p) is convex. □ 

Proof of Theorem I3.lt (i) We only need to consider the existence since the uniqueness is 
a consequence of the convexity of E 2 d{^/p) m Lemma 13.41 For existence, we may apply the same 
arguments in Theorem 12. 11 where instead, for sequence p n = (<I> n ) 2 , we have to show 

liminf H(p n ) > H(p°°), with p°° = |$°°| 2 . (3.17) 

n— yoo 

Denote 

<p n = (d n±n± - njA) [(-A)-VV], n = 0, 1, . . . , or n = 00. 



Using $ n -> $°° in L 2 (R 2 ) and <J>™ <J>°° in ff^R 2 ), then p™ -> p°° in LP(R 2 ) p > 1, and Lemma 
shows that -> in TU~ 1,P '(M 2 ) which is the dual space of with p' = p/(p - 1). Thus (|3~T7| 
is true and the existence of ground state follows. 

(ii) To prove the nonexistence results, we try to find the case where E 2 o doesn't have lower bound. 
For any $ e S 2 , denote p(x) = l^fx)! 2 and let 9 € R such that (cos#,sin#) = 1 (n t re 2 ) when 

n i + n i 7^ an d ^ = if rti = ri2 = 0. For any ex, £2 > 0, consider the following function 

$ £li£2 (x,y) = e 1 1 ^ 2 e 2 " 1 ^ 2 $(e^ 1 (a;cos0 + y sin9), e 2 x (—x sinO + ycos9)), (3.18) 
denoting p ei ,e 2 = |$ EljE2 | 2 , then 

P^(£,i, £2) = p(£i(Ci cos6» + £ 2 sin6>), e 2 (-£x sin6> + £ 2 cos0)). (3.19) 
By the Plancherel formula and changing variables, we get 

H(P^) = ^ J r2 ^ |P £1 , £2 | di 

= A/ ^ + ^^-^i 2 |^ (£1??1 , £27?2)d7? 

4tt 2 J R 2 1 7? I 

\P\ {VWn2)dr). 
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Let K = §^ , then the dominated convergence theorem implies 

H{Pe lt e 2 ) = < _„2 + 1 o(1) 2 (3.20) 

For fixed k > and letting ej — > + , we have J R2 ^(x)|$ £li£2 1 2 rfx = 0(1) and 

l|V# ei ,e,||l = -iR $ ll2 + 4ll»V*H2, ll*s liea ||1 = — ||*|lt (3-21) 

Thus under the condition (Bl"), i.c n% ^ and A > 0, choosing k large enough, we get 

= ^ + ^2+^2+ C 4 A ^ i+ 3 0(1) + 0(1), (3.22) 

where (fc = 1,2,3,4) are constants independent of k, e± and C4 > 0. Since 713 7^ 0, the last term 
is negative for n large, sending E\ — > + , one immediately finds that lim E 2 D(^e 1 .e 2 ) = —00, 

which justifies the nonexistence. Under the condition (B2"), i.e. n\ < ~ and A < 0, by choosing k 
small enough in (|3.20p , sending e\ to + , we will have the same results. Case (B3") will reduce to 
Theorem O □ 



3.2. Existence results for the Cauchy problem. Let us consider the Cauchy problem of 
equation (|1.15|) . noticing the nonlinearity (f>(d xl±Yl± — n|A)((— A) _1 / 2 |0| 2 ) is actually a derivative non- 
linearity, it would bring significant difficulty in analyzing the dynamic behavior. The common approach 
to solve the Schrodinger equation is trying to solve the corresponding integral equation by fixed point 
theorem. However, the loss of order 1 derivative due to the nonlocal term will cause trouble. This 
can be overcome by the smoothing effect of inhomogeneous problem iu t + Am = g(x, t), which pro- 
vides a gain of order 1 derivative [TUJ[T7]. When V 2 (x) = 0, i.c. without external trapping potential 
which corresponds to the free expansion of a dipolar BEC after turning off the confinement, the above 
approach can be extended straightforward. However, when V^(x) 7^ 0, i.e. with an external trap- 
ping potential, especially a confinement trapping potential with limixi^oo V^(x) = 00, the approach in 
[TUl [T7] has some difficulties. By configuring that (9 n±ni — n|A)((— A) _1 / 2 |</>| 2 ) is almost a first order 
derivative, we are able to establish the well-posedness of (| 1 . 1 5[) with a general external potential V% (x) 
by a different approach. 

The Cauchy problem of the Schrodinger equation with derivative nonlinearity has been investigated 
extensively in the literatures |161 118] . Here, we present an existence results in the energy space with 
the special structure of our nonlinearity, which will show that the approximation (|1.15j) of (|l.llj) is 
reasonable in suitable sense. We are interested in the case of A 7^ 0. 

Theorem 3.5. (Existence for Cauchy problem) Suppose the real potential V<2(x) satisfies \2. 7| ) 
and lini|x|->-oo U 2 (x) = 00, and initial value </>o(x) £ X2, either condition (B2) or (B3) in Theorem \3.1\ 
holds with constant Cb being replaced by C& / 1 1 ^"o 1 1 2 j then there exists a solution 4> G L°°([0, 00); X 2 ) H 
W '""([Qj 00); X%) for the Cauchy problem of \1.15\) . Moreover, there holds for L 2 norm and energy 
E2D i3.1\) conservation, i.e. 

IM-,t)||i«(R a ) = ||fo|U»(R«), E 2D {<f>(;t))<E 2D (fo), Vi>0. (3.23) 

Proof. We first consider the Cauchy problem for the following equation, 

0t^W) =H^ 5 + 5l (0 5 )+ 52 (/), xel 2 , <>0, (3.24) 
with the initial data 0*(x, 0) = 0o(x), /3 = tz^+B^- : ip s = Uj D * \<j) S \ 2 where Uj D is given in (HT2l 



2ir e 



as t/ 2D (x) = J r ^ +s2g2 ds (5 > 0), and 

1 3A 

H£ = -±A + V 2 (x), = /3 |/| V, ff2 (/) = -y/(9 n ^-7 ^ 2 A)/. (3.25) 
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Then our quasi-2D equation II (|1.15p can be written as 

^ = H^> + 5l (0) (3.26) 

where 

U<P) = ~ T <K9n±n ± ~ ^A)(-A)^/2(| (/) |2 ) _ (3 37) 

We denote the pairing of X 2 and its dual X 2 by (,}x 2 ,x% as 

(h, h)x 2 ,x* 2 = Re / /i(x)/ 2 (x) dx, (3.28) 
Js. 2 

where Re(/) denotes the real part of /. Using the results in Theorem 12.41 and [T3"] . we see that there 
exists a unique maximal solution <p s e C([— T^ in , T^ ax ], X 2 ) n C 1 ([— T^ in , T,* ax ], A"|). Here maximal 
means that if either t f T^ ax or 1 1 -T^ in , ||0' 5 (t)||x 2 -> oo. We want to show that as (5 -> 0+, (/i" 5 will 
converge to a solution of equation (|1.15l) . 

First, we show that T^ in = +00 and X^ lax = +00. The energy conservation for p.24[) is 

E s (t):= i||V^||! + l8b||/ll4+ / ^ 2 (x)|^| 2 dx + ^ ip (i)=£; a (0), i > 0, (3.29) 

where 

E S dip (t) = --r \<t> S \ 2 (d a±a± -n 2 A)/dx, *>0. (3.30) 

Similar computation as in Lemma 13.41 confirms that E^ ip > 0. Hence energy conservation will imply 
that U s {t)\\x 2 < 00 for all t, i.e. T n 5 lax - T 5 min = +00. 
We notice that 

X 2 ^ H 1 ^ L 2 ^ H- 1 ^ X;, (3.31) 

where is viewed as the dual of H . Consider a bounded time interval I = [— T, T], it follows from 
energy conservation that there exists a constant Ci(c/>q) > such that 

|| <P S || C([-T,T];X 2 ) < Cl(0o)- (3.32) 

Moreover, Lemma |2~21 and Remark |3. II would imply 

U s (d n±n± -r%£)<p% <C||^|| 9 .||V|^| 2 || P <C||^|| 9 -||0 5 ||2 P /(2- P )l|V^|| 2 , (3.33) 
for q,p £ (1, 2), i + i = i. Then we have 

||/||C1([-T,T];X 2 «) < C 2 (0o). (3-34) 

Thus, from (|3.32l) and (|3.34l) . there exists a sequence 5„ — > + (n = 1,2,...,) and a function <f> £ 
L ca ([-T,T};X 2 )^W 1 ' co {[-T,T];X*) [13] , such that 

4> 5 "{t) ^<p(t) in X 2 , for alU e [-T,T]. (3.35) 

For each t 6 [— T, T], due to the mass conservation of the equation p.24j) . we know ||<^ 5 "(i)||2 = ll^olh, 
by a similar proof in Theorem 12. 1[ the weak convergence of (i) in X 2 would imply that ^ (t) 
converges strongly in L 2 , which is a consequence of the fact that V 2 (x) is a confining potential. So, 
lim \\(/) Sn (t)\\ 2 = || (j){t) || 2 , and it turns out that [13] 

n—>oo 

in C([-T,T];L 2 (K 2 )). (3.36) 
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In view of ()3.35p . (|3.36|) and the Gagliardo-Nirenberg's inequality, we obtain 

in C([-T,T];L P (R 2 )), for aU p G [2, oo). (3.37) 

We now try to say that 4> actually solves equation (|1.15[) . For any function i/>( x ) G X2 and f(t) G 
C™([-T,T}), from equation we have 

' [<i0 s »,V)x 2 ,xj/(*) + (H^*» +9i(<P K ) +Sa(0 in ),^)x a ,xj/(t)] * = 0. (3.38) 

-T 

Recalling ^(u) - gi(w)| < C(|u| 2 + M 2 )|u - v\, (pH?) implies that 13 for all t G [-T, T] 

gi(4> Sn (t)) g x {(j>{t)), in L p (R 2 ) for some p G [l,oo), (3.39) 

W 6n (t)),i/>(t))x 2 ,x S -> (giW)),m)x 2 ,Xi- (3.40) 

For g2{<j) Sn ), consider y 5 "(x, t), noticing d a ip s ™ = T'q" ( | <^ 5 " | 2 ) (a = x,y) (defined in Lemma [3.3p . we 
have proven in Lemma 13.31 that T^ n is uniformly bounded from LP to LP and 

Tt(\m\ 2 )^R a (\m\ 2 ) = d a (-A)-^ 2 (\m 2 )^^(R 2 ), P 6(1,00), <5„^o+. (3.41) 

Rewriting 

^"(|^(t)| 2 ) - Tt{\^{t)\ 2 - |0(t)| 2 ) +Ti"(|0(t)| 2 ), (3.42) 
recalling the fact (|3.37p . we immediately have 

Tt(\<t> 5n (t)\ 2 ) -> ^(l^(t)l 2 ) in L P (K 2 ), for some p G (l,oc), (3.43) 

which is actually 

d a <p s "(t) -> d a ((-A) -1 / 2 !^)! 2 ) , inL p (K 2 ), for some p G (1, oo). (3.44) 
Hence, integration by parts, for a' = x, y, 

{(t> 5n {t)d aoi ,y 5 -(t),i>{t))x 2 ,x* 2 =Re / <f> Sn (t)d aa 'ip Sn {t)iP(t)dx 



= -Re / d a ^(t){d a ,<f> d »(t)i)(t)+<t> d -(t)d a/ ip(t))dx, 
Jr 2 

passing to the limit as rt — > oo, 

lim (<t> s "(t)d aa ><p s "{t)Mt))x 2 ,x S =-Re / i? Q (|0(i)| 2 )(a Q ^(t)^(t) + 0(i)9^?(t))dx 

= (#t)0«a' (-A)- 1/2 (|0W| 2 ), ^)>X 2 ,XJ , 

in view of (|3.44l) and (I3.35[) . we obtain 

lim (ff 2 (/"(t)),VW)x 3 ,X| = &(.<Kt))Mt))x 3 ,x S . (3.45) 

Combining the above results and (|3.40[) together, sending n — ► oo, dominated convergence theorem 
will yield 

T 

[(#, i>)x 2 ,x-J'{t) + (H^> + ffl (0) + .g 2 (0)»x 2 ,xj/(t)] df = 0, 

-T 
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which proves that 

t = H y > + < ?1 (0)+0 2 (0), in* 2 *, a.a. t G [— T, T], (3.46) 

with 0(f = 0) = 0o, and G L°°([-T, T]; A" 2 )nVt /1:00 ([— T, T]; X%). Moreover, by lower semi-continuity 
of X 2 norm, (|3.39|) and (|3.45l) . the energy i? 2 D (|3.ip satisfies 

E 2D {(t>(t)) <E 2D {fa). (3.47) 

It is easy to see that we can choose T = oo.D 

If the uniqueness of the L°°([—T, T];X 2 ) n W 1,0O ([— T, T]; A|) solution to the quasi-2D equation 
II (|1.15|) is known, we can prove that the solution constructed above in the Theorem 13.51 is actually 
C([-T, T];X 2 ) n C 1 ([-T, T]; X|) and conserves the energy. 

Next, we discuss possible finite time blow-up for the continuous solutions of the quasi-2D equation 
II (|1.15p . To this purpose, the following assumptions are introduced: 

(A) Assumption on the trap and coefficient of the cubic term, i.e. V 2 (x) satisfies 3y 2 (x) + x • 

VV^fx) > 0, P~ x +^ Xn 3 . > — ,,? b n -i , with 0o being the initial data of equation (11.151) : 

V ZTT € 1 1 I ! 2 

(B) Assumption on the trap and coefficient of the nonlocal term, i.e. V 2 (x) satisfies 2V% (x) + x • 
W 2 (x) > 0, A = or A > and n% > \. 

Theorem 3.6. (Finite time blow-up) For any initial data 0(x, t = 0) = 0o( x ) G A 2 with 
J R2 |x| 2 |0o( x )| 2 dx < oo, if conditions (Bl), (B2) and (B3) with constant Cb being replaced by Cb/||0o|| 2 
are not satisfied, let :— 0(x, t) be a C([0, T max ), X 2 ) solution of the problem \1.15\l with L 2 norm 
and energy conservation, then there exists finite time blow-up, i.e., T max < oo, if one of the following 
condition holds: 

(i) -E 2 _d(0o) < 0; an d either Assumption (A) or (B) holds; 

(ii) E 2 r){4>o) = and Im (J R2 0o(x) ( x ' V0o(x)) dx) < 0, and either Assumption (A ) or (B) holds; 
(Hi) E 2 d(M > 0, and Im (J* R2 O ( X ) (x • V0o( x )) dx) < -\j3E 2D ((/) )\\x(j) \\ 2 if Assumption (A) 

holds, or Im (J R2 0o( x ) ( x ■ V0o( x )) dx) < — ^/2i? 2 ,D(0o)||x0o|j 2 if Assumption (B) holds. 
Proof. Calculating derivatives of the variance defined in (|2.40p . for a = x,y, we have 



dt 

and 



— u a {t) = 2 Im I / ad a (j) dx ) , t > 0, (3.48) 



jL«7 a (t) = / [2|a Q 0| 2 + /?o|0| 4 + 3A|0| 2 aa Q (a„ ini - n 2 A)<^ - 2a|0| 2 d Q ^ 2 (x)] dx, (3.49) 
at J K 2 

where /3 = ^y^" 3 , (-A) 1 / 2 ^ = |0| 2 . Writing p = |0| 2 , (p = (3 n±ni - ti 2 AV and noticing that p 
is a real function, by the Plancherel formula, similar as Theorem 12. 5[ we get 



/ |0| 2 ( X -V^) dx = ~ [ |0| 2 ^dx. 

Jr 2 z Jm 2 



(3.50) 



Hence, summing p.49[) for a = x, y, and using the energy conservation, if Assumption (A) holds, we 
have 



^a v (t) = 2 J (JV0| 2 + /3 O |0| 4 - ^A|0| 2 (d n±n± - n 2 A) v - |0| 2 (x • V7 a (x))j dx 
= 6£ 2D (0)- / (|V0| 2 +/? o |0| 4 )dx-2 / |0| 2 (3V 2 (x)+x- W 2 (x)) dx 

JR 2 JR 2 

<6£ 2D (0(-,t)) <6£ 2D (0 O ), t>0. (3.51) 
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Thus, 

<rv(t) < 2,E 2D {(j) Q )t 2 + a' v {0)t + cry(O), t > 0, 



and the conclusion follows as in Theorem l2.5l If Assumption (B) holds, the energy contribution of the 
nonlocal part is non-positive and we have 



dt 2 



a v {t) = 2 J (jV0| 2 + /3oH 4 - ^AH 2 (d n±n± - n 2 3 A) ^ - |0| 2 (x • W 2 (x))J dx 
= 4E 2D {0)-^ [ \<f)\ 2 <fdx-2 [ |0| 2 (2F 2 (x)+x- W a (x)) rfx 

Z JR 2 JR 2 

<4E 2D (<f>(;t))<iE2D(<l>o), t>0, (3.52) 

and the conclusion follows in a similar way as the Assumption (A) case. □ 

4. Results for quasi-lD equation. In this section, we prove the existence and uniqueness of 
the ground state for quasi-lD equation (jl.171) and establish the well-posedness for dynamics. 

4.1. Existence and uniqueness of ground state. Associated to the quasi-lD equation (|1.17p . 
the energy is 



Bld($) 



i|3 2 <i>| 2 + Vl (z)\<P\ 2 + i/3 1D |$| 4 + 3A(1 J^W y 
2 2 16v27T£ 



dz, (4.1) 



where (3id = l3+2 ^^ 2 3n ^ anc j 

1 Z" 00 2 

<p(z) = d zz {Ul° *l*| 2 ), Ul D {z) = ^ e'^ds. (4.2) 

V 71 " 

Again, the ground state $ ff € Si of the equation (| 1 . 1 7[) is defined as the minimizer of the nonconvex 
minimization problem: 

Find $ g € Si, such that E 1D {<k g ) = min E 1D ($). (4.3) 
For the above ground state, we have the following results. 

Theorem 4.1. (Existence and uniqueness of ground state) Assume < V\{z) € Lj^ c (K) and 
lim Vi(z) — oo ; for any parameter j3, A and e, there exists a ground state <I> 5 € S\ of the quasi- 

\z\— >oo 

ID equation ^1.17}) - [1.18\) . and the positive ground state is unique under one of the following 

conditions: 

(CI) A(l - in 2 ) > and? - (1 - 3nf)A > 0; 
(C2) A(l - 3n 2 ) < and (i + f (1 - 3n§) > 0. 
Moreover, § g = e ie °\$> g \ for some constant 9q € R. 

To complete the proof, we first study the property of the convolution kernel U\ B (|1.18[) . 
Lemma 4.2. (Kernel Ul D il.l~S\) ) For any f(z) in the Schwartz space 5(R), we have 

- - — - \/2ft~(£) f°° e - e2s / 2 

JlD * fff\ — f(t\TnD(£\ — v *^-> I 

Hence 



Ul D * /(£) = muu>® = ^' I J^-/-^ ? e K. (4.4) 



|^(^*/)|| 2 <^2||/|| 2 . (4.5) 

V7T£ 
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Proof. For any f(z) G <S(R), rewrite the kernel as [TT] 



-I,:*,. , V?e f w 2 £ (x,y)w 2 (x',y') 



u e ( z ) = —7= i 1 : % ' * , ' =dxdydx'dy', 
applying Fourier transform to both sides and using the Plancherel formula as in Lemma l2.2[ we obtain 

Then a direct computation yields the conclusion. □ 

Lemma 4.3. For the energy E\r>{') in j4-l\ )> we have 

(i) For any <3? G S\, denote p(z) — |$(z)| 2 , then 

E 1D ($) > E 1D (\<Z>\) = E (^fp) , V*€5i, (4.7) 

so the ground state $> g of &4-l\ ) is of the form e l6o \^ > g \ for some constant 8o G R. 

(ii) Eid is bounded below. 

(Hi) If condition (CI) or (C2) in Theorem \4-.l\ holds, E\r>i^Jp) is strictly convex in p. 
Proof. Part (i) is similar to that in Lemma |2"T2"1 Part (ii) is well-known, once we notice the property 
of kernel U\ D (cf. Lemma \A. 2 1) and the Sobolev inequality in one dimension, 

Wfto < ll/lhll/'lb- (4.8) 

(hi) We come to the convexity of Eio(y^p). Following Lemma [2.31 we only need to consider the 
functional 



Hid(p) 



+ A(l-3n§)/2 2 , 3A(l-3n§) _ IJjXD 



4tt£ 2 16V2tt£ 2 



-pd zz (U,*p) 



dz. (4.9) 



Then under condition (CI) or (C2), using the Plancherel formula and Lemma [472J after similar com- 
putation as in Lemma 12. 2\ we would have H\d{p) > 0- For arbitrary Jp~x-,Jp~i G S\ and 6 G [0, 1], 
denote p = 9pi + (1 — 0)p2, then ^fp~ e g S\ and 

6H 1D { Pl ) + (1 - 9)H 1D (p 2 ) - H 1D { Po ) = 9(1 - 9)H 1D ( Pl - p 2 ) > 0, (4.10) 

which proves the convexity. □ 



Proof of Theorem I4.lt The uniqueness follows from the strict convexity in Lemma 14.31 The 
existence part is similar as Theorem 12.11 and we omit it here for brevity. □ 

4.2. Well-posedness for the Cauchy problem. Concerning the Cauchy problem, Lemma l4~2l 
shows that the nonlinearity in quasi- ID equation (jl.171) is almost a cubic nonlinearity, while the same 
property has been observed for quasi-2D equation I (|1.11[) - (|1.12|) . Hence similar results as Theorem 
12.41 can be obtained for equation (|1.17j) and we omit the proof here. 

Theorem 4.4. (Well-posedness for Cauchy problem) Suppose the real-valued trap potential satis- 
fies Vx{z) > for z G R and Vx(z) G C°°(R), D k Vx{z) G L°°(R) for all integers k>2, for any initial 
data <p(z,t = 0) = (f>o(z) G X\, there exists a unique solution <fi G C ([0,oo), X\) n C 1 ([0,oo), AT^) to 
the Cauchy problem of equation ji.iTj ). 

5. Convergence rate of dimension reduction. In this section, we discuss the dimension 
reduction of 3D GPPS to lower dimensions. Inspired by the previous work of Ben Abdallah et al. [8, 9 
for GPE without the dipolar term (i.e. A = 0) and 0(7] for Schrodinger-Poisson systems, we are going 
to find a limiting £-independent equation as e — » + . Thus in quasi-2D equation I (jl.lip . II (|1.15p and 
quasi-lD equation (|1.17p . we have to consider the coefficients to be 0(1). The existence of the global 
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solution for the full 3D system (|1.4j) - (|1.5p has been proven in [HI [2] when /3 > and A £ 
hence we would expect the limiting equation in lower dimensions valid in a similar regime. Thus in 
lower dimensions, we require that in the quasi-2D case, ft = 0(e), A = 0(e), and in the quasi-lD case, 
j3 = 0(e 2 ), A = 0(e 2 ), i.e. we are considering the weak interaction regime, then we would get an 
e-independent limiting equation. In this regime, we will see that GPPS will reduce to a regular GPE 
in lower dimensions. 

5.1. Reduction to 2D. We consider the weak interaction regime, i.e., f3 —> e/3, A — > eA. In Case 
I (|1.8[) . for full 3D GPPS (|1.4[) - (|1.5[) . introduce the re-scaling z — > ez, iji — > e" 1 / 2 -^ 5 which preserves 
the normalization, then 



idtip s (x, z, t) 



Hi 



^H z + (/3-A)|^| 2 -3eAa nen ^' 



(x,z) el 3 , t> o, 



where x = (a;, y) £ R 2 and 



H 



1 



-(d xx + d yy ) + V 2 (x,y), 



1 z 2 
H, = ——d zz H — — 

2 2 



(ni, n 2 , n 3 /e), <9 ne = n s • V, <9 ne „ e = d ne (d ne ), 

{-d xx - d yy - \d zz )ip £ = -m 2 , Um ^ £ (x)=0. 
e e |x|— >oo 



(5.1) 



(5.2) 
(5.3) 
(5.4) 



-z 2 /2 



It is well-known that H z has eigenvalues fik = k + 1/2 with corresponding eigenfunction Wk(z) (k 
0, 1, . . .), where {u>/c}fcL f° rms an orthornormal basis of L 2 (R) [THUS], specifically, Wq(z) — ^rpr* 
Following [5], it is convenient to consider the initial data concentrated on the ground mode of H z , i.e., 

-0 e (x, z,0) = <fio(x.)w (z), 0o e X 2 and ||0o|U 2 (R2) = 1- (5.5) 

In Case I (11.81) . when e — » + , quasi-2D equation I (jl.lip . II (|1.15[) will yield an e-independent 
equation in the weak interaction regime, 



ift^(x,,t) = H^+ ^ (1 / _ 3n ' )A |^| 2 0, x = (as,y) e R 2 , 



>2it 



(5.6) 



with initial condition 0(x, 0) = 0o( x )- We follow the ideas in [HI [8, 9 to show the convergence from 
the 3D GPPS to the 2D approximation. First, let us state the main result. 

Theorem 5.1. (Dimension reduction to 2D) Suppose V 2 satisfies condition — § < A < j3 

and /3 > 0, let ijf £ C([0, oo); X3) and <f> £ C([0, 00); X 2 ) be the unique solutions of equations h5.1}) - 
~5\) and i5.6\) , respectively, then for any T > 0, there exists Ct > such that 



ip e (x, z, t) — e t ~? r <f)(x,t)wo(z) 



<C T e, Wte[0,T] 



(5.7) 



Under the assumption, we have the global existence of i/j s [21 US] as well as <f) [SJ [13] . Define the 
projection operator onto the ground mode H z by 



Uip £ (x,z,t) = e 



(t> e {x,t)w (z), 



(x, z) £ 



t > 0, 



(f> e (x,t) =e lt *° t/e2 ip E (x,z,t)w a (z)dz, x = (x, y) £ R 2 , t>0. 



where 



Since the space (x, z) £ R 3 is anisotropic, we introduce the Lf space by the norm 
11/11 0,<z) := = II II/(-' z )IIl|IIl?, [l,oo]. 



(5.8) 
(5.9) 

(5.10) 
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The corresponding anisotropic Sobolev inequalities are available [8]. 

Lemma 5.2. (Uniform bound) Let tp £ and 4> be the solutions of \5.1\) and i5.6\) . respectively, <ff 
be defined in h5. 9\) . A e [— §,/3] and j3 > 0, we h 



ave 



V/ € L oo ((0,oo),iJ 1 (R 3 )), cj),(f> e e ^((O,! 1 ),^ 1 ^ 2 )), (5.11) 
with uniform bound in e. Moreover, for p £ [2, oo], 

\\r - n^ii 2 L2(R3) + \\e z (r - < ce\ \w - ww^ < c e , (5.12) 

with C depending on ||<^o||x 2 ; uniform in time t. 

Proof. From energy conservation of equation (|5.1[) , we know that the following energy is constant, 



E(t) := (H^ £ (t),^W) + ? (H^ e («),V £ (*)) + ^\m\h m + — II^V^WHl^.). 

where (•,•) denotes the standard L 2 inner product in K 3 . Using standard L p estimates for Poisson 
equation (|5.4I) . we have ||9„ e V<£ e (t)||i,2( R 3) < ■j||V' e (i)|lz,4( R 3), which implies 

+ ^-||9n £ V^(t)||| 2(R 3 ) >0, and £( ) = g+C , (5.13) 

oo 

where Co depends on ||</>o||x 2 - Writing ip e (t) = ^ (fik(x-,t)wk(z), and using the L 2 conservation 

fc=0 

oo 

||0fe(i)|! 2 2(- R 2' ) = 1, we can deduce from the energy conservation that 

k=0 



^ i 2 (R 2 ) 



1 1 

g + Co > (H^ e (*), V» e (*)) + ? (H,V £ (t), V> £ (*)) = CH^^(t),^(*)) + 3 5>*ll&(*)l 

1 OO 

= (H^ e (i),^(i)) + -5 - W>)ll&(*)ll W») + 



fc=0 



fe=i 



Hence, 



\\d x r\\h m + \\9 v r\\h m < mr(t),^(t)) < Co, (5.14) 

\\d z r\\l H m*) < (^r, r) < ^ + c o£ 2 , (5.15) 

\\r - n^|| 2 L2(R3) < — - IM>)Wk(t)\\l>W < 2C o^ 2 , (5.16) 



Mi - Mo k=1 



\\d z (r - n^ e )||i 2(R3) < ]T 77-^77^ - W>)ll^(*)lli'(R=«) < nC e 2 . (5.17) 

Estimate on \\ip € — ni/j e || ( p ,2) follows from Sobolev embedding. □ 

We also need the following Strichartz estimates for the unitary group e x , which is valid when 
V2 satisfies condition (|2.37[) [T5] . 

Definition 5.3. In two dimensions, let q' , r' be the conjugate index of q and r (I < q, r < oo) 
respectively, i.e. 1 = 1/q' + 1/q = X/r' + 1/r, we call the pair (q,r) admissible and (q' ,r') conjugate 
admissible if 

2 2(i-i), 2<r<oo. (5.18) 



q \2 r 
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The following estimates are established in |T3l [12j [26] . 

Lemma 5.4. (Strichartz's estimates) Let (q,r) be an admissible pair and (7, p) be a conjugate 
admissible pair, Jclfefl bounded interval satisfying 6 J, then we have 

(i) There exists a constant C depending on I and q such that 



< C(/,g)||v?|| L 2 (R2) . 



Li(I,L*-(M?)) 

(ii) If f € L^(I, L P (K 2 )), there exists a constant C depending on I, q and p, such that 



-i(t-s)H* 



f(s)ds 



< C(I,q, p)\\f\\ L -y(I <L r(Jg?)). 



(5.19) 



(5.20) 



Li(I,L r (R 2 )) 



Now, we are able to prove the theorem. 

Proof of Theorem I5.lt In view of Lemma r5.2[ we can derive 

\\i/) e - e _l ^"(/)Wo(z)||L2 (R 3) < \\ip £ - I1V' £ ||l2( R 3 ) + ||IT/; £ - e~ l ^~ </>u> (z)|| L 2 ( 

<C£+\\(f> s (t)-(t>(t)\\ LHA2) . 



(5.21) 



Hence, we need to estimate the difference between cjf and <j>. By (|5.ip and (|5.4p . we know </> e (x, t) in 
(|5.9|) solves the following equation 



id t (j) e = + (p - A + 3nl\)e ltM}t/E2 / \ip s \ 2 ip e w (z)dz + eg e , 

Jr 

g" = e^ 2 [ P s (<p*)rw (z)dz, 

JR 

where the differential operator P £ is defined as 

2 

Pe{v £ ) = -3A((n! - n 3 )9 xa; + (n 2 - ".3)^2,2, + 2n 1 n 2 d xy + -{n x n 3 d xz + n 2 n 3 d yz ))<p e . 

Denote x e (x,i) = (ff — cf), noticing that 1 1 tc?o 1 1 4 — 1/V27T, x e satisfies the following equation 

X e = H£ x e + f{ + /| + eg", x £ (* = 0) = 0, 
/3- A + 3n|A 



(5.22) 



/| = 08 - A + 3^A)e^*/ £2 / (|^ e |V - e-^*/ £2 |0 £ «, o | V™o) Wo(*)«k- 

JR V ' 

Applying Strichartz estimates on bounded interval [0,T] and recalling that (00, 2) is an admissible 
pair, we can obtain 

IIX £ ||l~([0,T];L 2 (R 2 )) < C [||/i ||lp*([o,T];Lp(R 2 )) + ll/f IIl->*([0,T];L-'(R 2 )) + £ ll3 e |lL9*([0,T];L9(R 2 ))] j 

where (p*,p), (7*, 7) and (q*,q) are some conjugate admissible pairs. By a similar argument in [5], 
we have the estimates for f[ and /f which comes from the cubic nonlinearity, i.e. for appropriate 
pe (1,2) and 7 e (1,2), 

ll/l ILp*([0,T]:Lp(R 2 )) < C|IX"IIlp*([0,T];L 2 (R 2 )): II /I II Lt* ([0,T];L^(R 2 )) < C £ - (5.23) 

The basic tools involved are the Holder's inequality, Sobolev inequalities and the estimates in Lemma 
15.11 and we omit the proof of this part here for brevity. So, 



IIx1l=°([0,T];L 2 (R 2 )) < C(IIX £ Hlp*([0,T];Lp(R 2 )) + £||g e ||L<!*([0,T];L9(R2)) + £ )- (5.24) 
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Next, we shall estimate g £ . Let ip\ and ip\ be the solution of the rescaled Poisson equation (|5.4p 
with \ip £ \ 2 replaced by |H0 e | 2 and \ip e \ 2 — |n?/> e | 2 , respectively, then ip £ = <p\ + ipf, and we can rewrite 



<f = Jf + Jf + J: 



3 I 



(5.25) 



where 



Jf= / P^lWwidz, Jf =e 



•fa* 



P e (p £ )(ijj s - n^ £ )w dz, Jf = e^ 2 " / F £ (^)nv/worfz. 



For J £ , this one reduces to the quasi-2D equation I (II. lip , where we have that 
Jf = -3X(d n±n± - nfA)^^, and tf D = U 2D * \cj> £ \ 2 , 



(5.26) 



with U 2D given in (|1.12p . In view of the property of U 2D in Lemma 12.21 and Remark 13.11 recalling 
(ff 6 L°° ([0, oo); H 1 (M 2 )) , using Holder's inequality and Sobolev inequality, we obtain 



\J!\\ P < II^!d)IUI0 e IL < c||v|<f l 2 IUII<f IU < c, 



(5.27) 



where 1 < p < pi < 2, - = — + — . 

For Jf, applying Minkowski inequality, Holder's inequality and Sobolev inequality, as well as 
estimates for Poisson equation, noticing -0 £ S L°°([0, oo); i? 1 (K 3 )) and Lemma [5T2j we estimate 

\il> £ \ 2 

\\Ji\\ P < \\P s (v e )(r - n^>oll(i,p) < c\\p e (<p £ )\\ LP * {WL3) \\r - nr\\(oo,2) < CeII^-H^.^,) < c, 

where p* = 2p/(2-p) < 3. 

For Jf , similar as J*, J 2 , we have 

II -flip < \\PMW«>0\kl,p) < C\\Pe(<p e 2 )\\Lvi(R3)\\<t> e \\L*HW) < f II " IW 2 || LP1 ( R 3) 

c 

< —\\ip e - n^|| i 2 (IR 3 ) (||^ e || iJ3 3(R3) + \\n^ £ \\ L p 3(&3) ) < c, 

where p 3 = 2p\j [2—p\) < 6. Hence, by choosing p = 6/5, and p\ =4/3 would satisfy all the conditions 
for J| (k= 1,2,3), where we shall derive that uniformly in t, 



||s £ IUp(» 2 ) < II^iIUp(r 2 ) + ll^f IUp(r 2 ) + II^IIIlp(r 2 ) < c. 

Then choose q = p in (|5.24p , we have 

IIX e ||L~([0,T];L2(K2)) < C [||x £ ||lp*([0,T];L 2 (R 2 )) + e ] • 

Applying the results for all t € [0, T], we find 



llx e (*)ll£ <c 



||x e («)||£ ds + e"* 



te [0,T], 



(5.28) 



(5.29) 



(5.30) 



and Gronwal's inequality will give that |jx e (£)l|2 _ Ge for all t £ [0, T]. Combining with (|5.2ip . we can 
draw the desired conclusion. □ 

5.2. Reduction to ID. In this case, we again consider the weak interaction regime j3 —>■ e 2 /3, 
A -> e 2 A. In Case II fO]) . for the full 3D GPPS (jO)) -([l~5 |) . introducing the re-scaling x4 ex, y -> ey, 

■0 — > e _1 i/) £ which preserves the normalization, then 



id t ip £ {yi,z,t) 



r H x + (/3-A)|^| 2 - 3e\d Refle v 



*p e , (x,z)G 



(5.31) 
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where x = (x,y) € R 2 and 
n e = 

, 1 



1 



d zz + V 1 {z), 



1 



(d xx + d yy + x 2 +y 2 ), 



(n 1 /e,n 2 /s,n 3 ), d he 



1 



^dyy-d ZZ )ip E 



n e • V, dn €ilc = 9fi £ (9 fie ), 
1 

72 



£|2 



lim ip £ (x) 

|x|— >-oo 



(5.32) 
(5.33) 
(5.34) 



Note that the ground state mode of H x would be given by wo(x)wo(y) with eigenvalue 1, and the 
initial data is then assumed to be 



^ e (x, 2,0) = <j) (z)w (x)wo(y), (f>o e Xi and ||</>o||l 2 ( 



1. 



(5.35) 



In Case II (|1.9[> in the introduction, when e — > + , the quasi-lD equation (j!.17|) will lead to an 
e-independent equation in the weak interaction regime, 



id t (f>(z,t) 



^ + |A£-3n|) 
2^ ^ q 



z e 



t > 0, 



(5.36) 



with the initial condition 4>(z, 0) = 4>o{z). 

Following the steps in the last subsection, we can prove the following results. 

Theorem 5.5. (Dimension reduction to ID) Suppose the real-valued trap potential satisfies 
Vx(z) > for z <E R and V%(z) G C°°(R) 7 D^z) G L°°(R) for all k > 2. Assume -f < A < 
and /3 > 0, and let -0 £ € C([0, oo); X 3 ) anrf G C([0,oo);Xi) be the unique solutions of the equations 
i5.31\) - l[5. 35\) and i5.36}) , respectively, then for any T > 0, there exists C't > such that 



i(j £ (x,z,t) - e tt/e2 (t){z,t)wo(x)wo(y) 



L 2 (R3) 



<C T e, VtG[0,T]. 



(5.37) 



6. Conclusion. We have analyzed the effective lower dimensional models for three dimensional 
Gross-Pitaevskii-Poisson system (GPPS) describing dipolar Bose-Einstein condensates (BEC) in anisotropic 
confinement. The quasi-2D approximate equations I (II. lip and II (I1.15[) are introduced in the case 
that the trap is strongly confined in the vertical z-direction, and the quasi-lD approximate equation 
(I1.17P is presented in the case that the trap is strongly confined in the x-, y- directions. Properties of 
ground states for all equations, such as existence and uniqueness as well as non-uniqueness results were 
studied. Well-posedness of the Cauchy problems for both equations and possible finite time blow-up 
in 2D case were discussed. Finally, we rigorously proved the linear convergence rate of the dimension 
reduction from 3D GPPS to its quasi-2D and quasi-lD approximations in the weak interaction regime, 
i.e. P = A = 0(e 3 ~ d ), in lower d (d = 1,2) dimensions. In such situation, all the nonlocal terms in 
the effective equation (ll.llj) . (|1.15p and (|1.17l) vanish, resulting in a regular GPE in lower dimensions. 
We remark that the results in the paper hold true for a larger class of confinements rather that the 
harmonic ones. In fact, effective two-dimensional models have been derived and analyzed recently for 
a multilayer stack of dipolar BEC formed by a strong lattice potential 22 . 
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